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We investigate low-temperature transport properties of thin TiN superconducting films in the 
vicinity of the disorder-driven superconductor-insulator transition. In a zero magnetic field, we 
find an extremely sharp separation between superconducting and insulating phases, evidencing a 
direct superconductor-insulator transition without an intermediate metallic phase. At moderate 
temperatures, in the insulating films we reveal thermally activated conductivity with the magnetic 
field-dependent activation energy. At very low temperatures, we observe a zero-conductivity state, 
which is destroyed at some depinning threshold voltage Vt- These findings indicate formation of 
a distinct collective state of the localized Cooper pairs in the critical region at both sides of the 
transition. 



An early suggestion that tuning disorder strength can 
cause a direct superconductor-insulator transition (SIT) 
in two-dimensional systems jl| triggered explosive activ- 
ity in experimental studies of superconductor films Q. 
Experimentally, the SIT can be induced by decreasing 
the film thickness 3 and/or, close to the critical thick- 
ness, also by the magnetic field These scenarios 
are commonly referred to as disorder-driven SIT (D-SIT) 
and magnetic- field driven SIT. Recent studies on the B- 
induced insulator revealed several striking features: a 
magnetic-field-dependent thermally activated behavior of 
the conductivity and a threshold response to the dc 
voltage , indicating the possible formation of a distinct 
collective insulating state. Importantly, these findings re- 
fer to the films belonging to the superconducting side of 
the D-SIT. This rises the question of whether the above 
findings are specific only to the superconducting side of 
the D-SIT or a characteristic feature of the whole critical 
region including both the insulating and superconducting 
sides of the D-SIT. 

In this Letter we focus on the insulating side of 
the disorder-driven superconductor-insulator transition 
in TiN films. The transition itself turns out to be ex- 
ceptionally sharp. At zero and low magnetic fields we 
find thermally activated behavior of the conductivity. A 
positive magnetoresistance and a distinct threshold be- 
havior in the low-temperature I-V characteristics persist 
on the insulating side of the D-SIT. Our results clearly 
indicate that, in the vicinity of the D-SIT, the response 
to applied magnetic and/or electric fields, is the same 
irrespective of whether the underlying ground state is 
superconducting or insulating. 

The 5-nm thick TiN films were grown by atomic layer 
chemical vapor deposition onto a Si/Si02 substrate. The 
samples for transport measurements were patterned into 
Hall bridges using conventional UV lithography and sub- 



sequent plasma etching. To increase sheet resistances 
(i?n) without introducing structural changes, the films 
were thinned by an additional soft plasma etching. Elec- 
tron transmission micrographs and diffraction patterns 
revealed a polycrystalline structure in both initial and 
etched films, the interfaces separating densely-packed 
crystallites being 1-2 atomic layers thick. As we found 
before in such samples kpi — 1, where kp is the 
Fermi vector, and i ~ 0.3 nm is the mean free path. 
This short mean free path can be related to the enhanced 
CI content (up to 3%), characteristic of films grown by 
the above method [s!]. Four-probe resistance measure- 
ments were carried out by the standard low frequency 
(0.4-2 Hz) ac lock-in current source technique with the 
ac current 0.01-1 nA. In cases where the resistances were 
too high to employ the current source four-probe lock-in 
measurements, the two-probe voltage source technique 
with the ac voltage 10-30 fiV was used. At resistances 
R > 1 Mil, this allowed us to keep the power dissipation 
below 10^^^ W, thus ensuring linear regime and exclud- 
ing overheating. Two-probe differential conductance vs 
dc voltage measurements were done by means of the low 
frequency ac lock-in technique combined with the dc volt- 
age excitations. Magnetic fields up to 16 T were applied 
perpendicular to the film surface. 

We start with the zero magnetic-field results. Fig- 
ure [1]^ a) shows the temperature dependence of logi?n for 
seven different films. An increase in disorder results in 
a growth of i?n and reduces the superconductor critical 
temperature monotonically. We did not detect any sign 
of the reentrant behavior or a kink in i?(T), which are 
characteristic to the granular films and/or the films con- 
taining large scale inhomogeneities It indicates that 
our films do not have a granular structure but are rather 
homogeneously disordered. To characterize the behav- 
ior on the nonsuperconducting side, we replot the data 
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FIG. 1: (color online) . Temperature dependences of iZg taken 
at zero magnetic field for the samples near the localization 
threshold, (a) log Ra versus T. Inset: some part of the i?n 
data in a linear scale, (b) log versus 1/T for samples II, 12, 
and 13. Dashed lines represent Eq. (1) and fit perfectly at low 
temperatures. All curves saturate at the same 7?n ~ 20kQ 
at high temperatures, (c) Rn versus dashed lines are 

given by Rq = Ri exp{Ti/Ty/^ which (with Ri ~ 6kn) well 
fit the data at high temperatures. Vertical strokes mark To, 
determined by the fit to the Arrhenius formula of Eq. (1). 



for R\j versus 1/T in Fig.[IJb). At low temperatures we 
observe an Arrhenius behavior of the resistance, demon- 
strating that these samples are indeed insulators. The 
dashed lines correspond to 



R^Roexp{To/T). 



(1) 



The activation temperatures Tq in the three samples mea- 
sured are To = 0.25K (II), 0.38K (12), and 0.61K (13) 
(the growth in To corresponds to increasing disorder). 
The prefactor Rq determined from the extrapolation of 
the dashed lines in Fig.[lja) towards 1/T = is almost 
the same (w 20'k^l) for all samples. Figured] shows that 
our TiN films demonstrate an abrupt switch between 
the superconducting and insulating phases: indeed, the 
"last" superconducting and the "first" insulating films 
are practically indistinguishable by their resistances at 
temperatures higher than 1 K, for instance, at T = 10 K 
Ra = 8.74 kO (SI) and i?n = 9.16 kO (12). However, at 
lower temperatures they choose unequivocally between 
either the superconducting or insulating ground states. 

We point out that the separatrix between the sets of 
"superconducting" and "insulating" R{T) curves is not 
simply a horizontal line [R{T) = const] . To demonstrate 
that, we replot the data in the linear scale in the inset to 
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FIG. 2: (color online). Magnetoresistance isotherms for su- 
perconducting (SI, S2) and insulating samples (II, 13) at sim- 
ilar temperatures. All curves converge above 2 T. 



Fig. [Ha) and note that R{T) dependences of supercon- 
ducting samples are nonmonotonic. They exhibit an "in- 
sulating trend," i.e. an upward turn of the resistance pre- 
ceding its eventual drop to zero at low temperatures. Fur- 
ther insight into the evolution of TiN films across the D- 
SIT is drawn from the logi?n against 1/T^/^ plots shown 
in Fig. [Ijc). At T > To, the resistances of the insulat- 
ing samples compare favorably with the Efros-Shklovskii 

[HIIH. The tem- 



(ES) formula, R = i?i exp(Ti/T)i/2, 
peratures Ti in the three samples shown in Fig. [TJc) are 
Ti = 1.75 K (SI), 1.80 K (12), and 2.53 K (13). The pref- 
actor Ri is again nearly the same for all samples, but in 
this case it is close to the quantum resistance for pairs, 
/i/(2e)^ = 6.45 kil. At lower temperatures R{T) deviates 
from the ES behavior, which in the insulating samples 
turns into the Arrhenius law below To, while the super- 
conducting samples fall into a superconducting state. 

Turning to the magnetoresistance data shown in Fig. [21 
we see that in all samples, including the insulating films, 
R{B) varies nonmonotonically with B. It exhibits a posi- 
tive magnetoresistance (PMR) at low fields, then reaches 
a maximum, followed first by a rapid drop, and eventually 
saturates at higher magnetic fields [l2| , where the differ- 
ence between insulating and superconducting samples is 
suppressed and all curves converge. Since PMR in super- 
conducting films appears because of the suppression of 
superconducting phase coherence by the magnetic field, 
one can conjecture that this phase coherence persists also 
in insulating films. At low temperature the ratio of the 
magnitude of the resistance at maximum to its value at 
zero magnetic field increases. Figure [Hj presents R{B) of 
the sample II at low temperatures and shows that the re- 
sistance at low fields again follows the Arrhenius behavior 
with B-dependent Tq [Eq.(l)]. The activation tempera- 
ture Tq{B) qualitatively follows the magnetoresistance. 

We now discuss the most intriguing feature, a de- 
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FIG. 3: (a) Sheet resistance of sample II as a function of the 
magnetic field at some temperatures listed, (b) R versus l/T 
at B = (open circles), 0.2 (triangles), 0.3 (filled circles), and 
0.5 T (squares). The dashed lines are given by Eq. (1). (c) To 
(left axis), calculated from fits to Eq. (1), and the threshold 
voltage Vt (right axis) as a function of B. 
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FIG. 4: Differential conductance vs dc voltage for sample II. 

(a) d7/dV(Vdc) at B = 0.9 T at two temperatures T = 0.02 
and 0.07 K. An arrow shows the direction of the voltage sweep. 

(b) dlldV{Vdc) at B = 1.25, 1.65, 2.7, 6.2, and 16T at T = 
0.02 K. (c) d//d1/(Vdc) at B = 0, 0.1, and 0.2 T at T = 0.02 K. 



pinning transition observed at low temperatures in the 
insulating films. Figure |4] reveals an abrupt onset of 
finite conductivity as the bias voltage Vdc exceeds a 
threshold voltage Vt'- panel (a) presents a differential 
conductance, dl/dV versus Vdc measured on sample II 
at B = 0.9 T, using the two-probe technique with the 



contact separation 1.5 mm. Two traces represent data 
taken at T = 70 mK and at the lowest temperature, 
T = 20 mK, achieved in the experiment. The dl/dV 
curve at T = 70 mK is typical for an insulator with acti- 
vated conductance, showing a gradual increase of dl/dV 
with Vdc, and is symmetric with respect to the direction 
of the Vdc sweep. However, as temperature is decreased 
to r = 20 mK the voltage response changes dramatically. 
Under low Vdc, both, the current and dl/dV are immea- 
surably small. As soon as Vdc reaches some well-defined 
threshold value Vt, dl /dV abruptly jumps up by several 
orders of magnitude. The threshold behavior is accom- 
panied by a hysteresis. A sharp conductance jump is 
observed up to i? 2 T. A non-Ohmic behavior remains 
even at i? = 16 T. The threshold voltage changes non- 
monotonically upon magnetic field [see Vt{B) along with 
To{B) in Fig.[3l^c)]. Note the large magnitude of thresh- 
old voltage to activation energy ratio: eVT/kBTo « 220 
at B = 0.7 T. 

Notable is also the nonmonotonic behavior of 
dI/dV{Vdc) at B = 0, which is displayed in Fig. Hljc). 
Similar to the depinning-like behavior at finite magnetic 
field in Figs. [Ha) and (b), we find a steep initial in- 
crease of dl /dV{V) with a maximum around Vdc ~ 1 mV 
(marked by an arrow), followed by a gradual decrease. 
At the peak, dl jdV is about twice larger than at Vdc = 
4 mV. This feature cannot be explained by electron heat- 
ing and vanishes already at B = 0.025 T. Such a suppres- 
sion of the conductivity by the bias current or voltage is 
typical for a superconductor in a fluctuation regime. The 
observation a similar "superconducting trend" in the in- 
sulating regime indicates that superconducting correla- 
tions survive also on the insulating side of the D-SIT. 
This implies the presence of the localized Cooper pairs 
and local superconducting phase coherence. 

To summarize the essential findings in the D-SIT crit- 
ical region, we note that the search for a disorder-driven 
superconductor-insulator transition has included many 
materials^.g., Bi Q, MoSi [H, Ta [H, InO^^ SEIlJ, 
and Be [l7|. However, the immediate onset of exponen- 
tial temperature dependence of the resistance, which con- 
clusively evidences the direct transition into an insulator, 
was found so far only in InO^; [sl. 11 ^ . [Til and Be [13] films. 
In Bi, MoSi, and Ta compounds [sl. llSl. \v\ a weak loga- 
rithmic temperature dependence of the resistance was ob- 
served on the nonsuperconducting side, either because of 
a possible intermediate metallic phase or since the tran- 
sition to the superconducting state occurs there at much 
lower temperatures. Our data on TiN unambiguously 
show a direct D-SIT with the nonhorizontal R(T) sepa- 
ratrix between the insulating and superconducting sides. 
An insulating trend (i.e. the upturn of the separatrix) 
can also be seen in the data on InOa; [l3| and Be [Tsj 
films. This implies that D-SIT cannot be described by a 
single-parameter scaling with the universal resistance at 
the transition. 
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The next important feature is that in all three mate- 
rials exhibiting D-SIT, TiN, InO^, and Be, 0, [H [B, 
[igj . the low-temperature activated behavior, R{T) oc 
exp(ro/r) transforms into a variable range hopping upon 
increasing temperature, contrary to common wisdom ex- 
pectations. The fact that the high-temperature behavior 
of the last superconducting sample, SI, is close to that 
of 12 [R(T)s are nearly indistinguishable at high tem- 
peratures; see Fig. 1] evidences the insulating features of 
superconducting samples in the critical region. On the 
other hand, the similarity in R(T) between superconduct- 
ing and insulating samples in the ES regime indicates the 
presence of Cooper pairs in the insulating samples as well. 
Further, all these materials exhibit positive magnetoresis- 
tance on the insulating side of the D-SIT (InOa^ and 
Be films [l^Hlj). And, finally, nonmonotonic dependen- 
cies of Tq{B), extracted from the Arrhcnius behavior of 
R{T), and those of the threshold voltage, Vt{B), which 
we find at the insulating side of D-SIT, have also been 
observed in the magnetic-field-induced insulatingphase 
in samples, which are superconductors at zero B [1, 0]. 

We would like to emphasize that the investigated films 
are homogeneously disordered and do not possess struc- 
tural granularity. A positive magnetoresistance and 
threshold behavior, recently reported by [2^ . were ob- 
served in structurally granular systems at temperatures 
where granules are superconducting, while disappearing 
with the breakdown of superconductivity in granules. 

From all the above we conclude that in the critical re- 
gion of the transition a peculiar highly inhomogeneous 
insulating phase with superconducting correlations, a 
Cooper-pair insulator, forms. In other words, the films 
in the critical region may be regarded as possessing a 
self-induced granularity due to strong mesoscopic fluctua- 
tions in disordered superconducting thin films ^23j giving 
rise to formation of superconducting droplets (islands) 
immersed into an insulating matrix. Such a state can 
be viewed as a network of superconducting islands cou- 
pled by weak links (Josephson junctions arrays). This 
is strongly supported by experiments of (23 |. where the 
voltage threshold similar to ours was observed on the long 
chains of SQUIDs. 

The Cooper-pair insulator establishes as a result of 
the mutual Josephson phase locking and exhibits generic 
collective behavior [25j. At moderate temperatures, it 
shows thermally activated conductivity governed by the 
large collective Coulomb blockade gap for a Cooper-pair 
propagation. The observed voltage threshold behavior 
suggests that the Cooper-pair insulator falls, at very low 
temperatures, into a distinct zero- conductivity state. A 
nonmonotonic magnetic-field dependence of Tq and Vr 
(Fig. [3]) naturally results from the magnetic-field mod- 
ulation of the effective Josephson energy [2^ [2^. The 
appealing task now is the construction of the phase di- 



agram of the zero-conductivity state and revealing the 
mechanisms of depinning and temperature crossovers. 

We thank D. Weiss and W. Wegscheider for access 
to their high magnetic field system, V. F. Gantmakher, 
M. Feigel'man, A. Finkelstein, and M. Fistul for use- 
ful discussions. This research is supported by the Pro- 
gram "Quantum macrophysics" of the Russian Academy 
of Sciences, the Russian Foundation for Basic Research 
(Grants No. 06-02-16704 and No. 07-02-00310), the 
U.S. Department of Energy Ofhce of Science under the 
Contract No. DE-AC02-06CH11357, and the Deutsche 
Forschungsgemeinschaft within the GRK 638. 



[1] M.P.A. Fisher and D.H. Lee, Phys. Rev. B 39, 2756 
(1989). 

[2] For a review, see, e.g., A. Goldman and N. Markovic, 

Physics Today 51, No. 11, 39 (1998). 
[3] D.B. Haviland et al, Phys. Rev. Lett. 62, 2180 (1989); 

Y. Liu et al, Phys. Rev. B 47, 5931 (1993). 
[4] A. F. Hebard and M. A. Paalanen, Phys. Rev. Lett. 65, 

927 (1990). 

[5] G. Sambandamurthy, L. W. Engel, A. Johansson, and 

D. Shahar, Phys. Rev. Lett. 92, 107005 (2004). 

[6] G. Sambandamurthy, L. W. Engel, A. Johansson, 

E. Peled, D. Shahar, Phys. Rev. Lett. 94, 017003 (2005). 
[7] T.I. Baturina et al, Pis'ma v ZhETF 79, 416 (2004) 

[JETP Lett. 79, 337 (2004)]. 
[8] A. Satta et al. Mat. Res. Soc. Symp. Proc. 612, D6.5 
(2000). 

[9] see e.g. H. M. Jaeger et al, Phys. Rev. B 40, 182 (1989); 

A. Frydman et al, Phys. Rev. B 66, 052509 (2002). 
[10] B. I. Shklovskii and A. L. Efros, Electronic Properties of 

Doped Semiconductors (Springer- Verlag, Berlin, 1984). 
[11] A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1975). 
[12] T.I. Baturina et al, Phys. Rev. Lett. 98, 127003 (2007). 
[13] S. Okuma et al, Phys. Rev. B 58, 2816 (1998). 
[14] Y. Qin et al, Phys. Rev. B 73, 100505(R) (2006). 
[15] D. Shahar and Z. Ovadyahu, Phys. Rev. B 46, 10917 (1992). 
[16] V.F. Gantmakher et al, Zh. Eksp. Teor. Fiz. 109, 1765 

(1996) [JETP 82, 951 (1996)]. 
[17] E. Bielejec et al, Phys. Rev. Lett. 87, 036801 (2001). 
[18] W. Wu and E. Bielejec, cond-mat/051121. 
[19] D. Kowal and Z. Ovadyahu, Solid State Comm. 90, 783 

(1994). 

[20] E. Bielejec et al, Phys. Rev. B 63, 100502(R) (2001). 
[21] V. Yu. Butko and P. W. Adams, Nature 409, 161 (2001). 
[22] C. Christiansen, L. M. Hernandez, A. M. Goldman, Phys. 

Rev. Lett. 88, 037004 (2002). K. H. Sarwa B. Tan, 

K.A. Parendo, A.M. Goldman, cond-mat/0704.0765. 
[23] M. A. Skvortsov and M. V. Feigel'man, Phys. Rev. Lett. 

95, 057002 (2005). 
[24] E. Chow, P. Delsing, and D. B. Haviland, Phys. Rev. 

Lett. 81, 204 (1998). 
[25] M. V. Fistul, V. M. Vinokur, and T. I. Baturina, arXiv: 

0708.2334. 

[26] I. S. Beloborodov et al, Phys. Rev. B 74, 014502 (2006). 



